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Abstract 

Key words: point vortices, Riemann surfaces, Green functions. 
AMS MSC(2000) ; 76B47, 31C12, 37J05 . 

We consider N point vortices Sj of strengths Kj moving on a closed (compact, boundaryless, orientable) 
surface S with riemannian metric g. As far as we know, only the sphere or surfaces of revolution, the latter 
qualitatively, have been treated in the available literature. The aim of this note is to present an intrinsic 
geometric formulation for the general case. Since the pioneer works of Bogomolov [5] and Kimura/Okamoto 
[36] on the sphere S" 2 , it is known that stream function produced by a unit point vortex at s Q G S on a 
background uniform counter vorticity field is given by Green's function G g (s,s ) of the Laplace-Beltrami 
operator A g — divg o grad ff . It behaves as logd(s,s )/27r near s . Desingularizing the stream function 
G g is therefore mainstream mathematics (no pun intended). An energy core argument shows that a single 
vortex s drifts according to the Hamiltonian system (flg,R g ), where fl g is the area form of g and R g 
is Robins's function R g (s ) = lim s ^ So G s (s,s ) — logd(s, s )/2ir . If S has genus zero, it it known that 
R g = A~ l K + trace A^ 1 /A(S). The collective Hamiltonian keeps the same form of C.C.Lin's now classical 
paper [35] j with 

N N 
^collective = ^ KjftjSj) , H = ^ KiKjGgis^Sj) +^-KjRg(s e ) . 

3=1 l<i<j<N 1=1 

For Jordan domains D C S the structure is the same, using the appropriate hydrodynamical Green function. 
The extension for vortices with mass is also immediate. Under conformal changes of metrics g = h 2 g the 
symplectic form changes accordingly to fig = h 2 tt g while the new Hamiltonian is given by 

1 N K N N 

77 i=i a \") i=i i=i 

The presence of the total vorticity in the conformal change formula reflects the fact that when the sum of the 
vorticities vanish, the collective vortex stream function ip(s; s\, sjv) is independent of the conformal metric 
g = exp(2<p)g. In this case the ^-regularization of ip at any of the vortices simplifies: one just subtracts off 
^-<t>{_Sj) from the ^-regularized stream function at Sj, In particular, when S is conformal to the standard 
sphere, making an artificial puncture at any given s* £ S allows to easily write vortex motions for any 
metric on the sphere. We give a simple proof of Kimura's conjecture that a dipole describes geodesic motion. 
Searching for integrable vortex pairs systems on Liouville surfaces is in order. The vortex pair system on a 
triaxial ellipsoid extends Jacobi's geodesies. Is it Arnold-Liouville integrable? Not in our wildest dreams is 
another possibility: that quantizing a vortex system could relate with a million dollars worth question, but 
we took courage - nerve is more like it - to also present it. 
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1 Introduction 

The purpose if this note is to formulate vortex motions in intrinsic fashion. The (now classical) results of C.C.Lin 
([45], 1941) are generalized and reinterpreted. Complex coordinates z £ C will not be used (well, almost). The 
key observation is that the core energy desingularization procedure ([23, 1997) geometrizes in a nice way, matching 
perfectly with Green's function behavior. 

1.1 Vortex history on a capsule 

Helmholtz' "Wirbel" paper, short for ( Uber integrate der hydrodynamischen gleichungen welche den Wirbelbewegungen 
entsprechen, [31], 1858) came out a century+ after Euler's Principes generaux du mouvement des fluides ([25], 1757), 
where the mathematical definition of vorticity first appearecQ. Eighteen years later Kirchhoff presented the equations 
for point vortices in the plane (Vorlesungen, [37] . chapter XX, p. 259, eq.(14)): 

dxi 8P dyi dP D V- 1 
m ^ = dy~ ' m ^ = - ' P = ^ m * m > 

(nowadays vortex strengths are denoted by k or T/2n instead of m, 1/ir is omitted, and H replaces P). 

This is a Hamiltonian system (Q,P), where the symplectic form is the combination of planar areas weighted by 
the vorticities, Qsystem = ^2 Kjdxj A dyj . Only in 1941 the correct Hamiltonian was found by C.C.Lin, [45] . for the 
case of arbitrary planar domains D C C. See [49] and [7] for treatises and [T], [2], [6j for recent reviews. 

1.2 Vortices on curved surfaces 

Vortex dynamics on the sphere and hyperbolic plane started with Bogomolov [8] and Kimura and Okamoto [35] . 
Basically, since the late 1970's up to now only the constant curvature case has been considered. Follows a very 
incomplete list (we apologize for omissions): Kimura and Okamoto ([35], 1987), Pekarski and Marsden ([52], 1998), 
Kidambi and Newton ([33], 1998, [34], 2000), Lim et al. ([44] 2001), Cabral and Boatto ([H], 2003), Borisov and al. 
([ID]. [IT] . 1998, [12], 2000), Laurent-Polz 2005), Tronin ([61], 2006). The case of surfaces of revolution has just 

been touched recently, in qualitative terms ([57], 2002) or perturbatively ([16]. 2007). 

How about the general case? Vortex equations in the case of domains conformal to a planar domain were obtained 
by Hally [21]. The key aspect is how to desingularize the velocity field at a vortex. Using isothermic coordinates, and 
after some heuristics (complicated for mathematicians) he found additional "self-terms" attributed to the variations in 
the surface curvature near each vortex. For a surface homeomorphic to the z-plane C, with metric g — h 2 (z,~z)dz®d~z, 
Hally's equations are 

h 2 {z h Zj)% =i V] — + iKj-^-log(h(zj^j)) , (1 < j < N) . (1.1) 

Z£ Z -j CJZ<i 

A single unit vortex obeys the "self-motion equation" 

j = th- 2 (z,z)^-log(h(z,z)) . (1.2) 

We provide corrections to [21] for the case of nonzero total vorticity. Our results are valid for closed surfaces of 
arbitrary genus and Jordan domains on them. 

1 Circa 340 B.D. Aristoteles tried to explain typhoons in his Metereologica. The importance of vorticity in microfluid 
motion was clearly recognized by Leonardo da Vinci. Decartes envisaged vortices as powerhouses for planetary 
movement. W. and J.J. Thomson used vortices for an atomic theory, discarded after Rutherford's 1912 experiments. 
Some say string theory is a souped up version of the primeval vortex soup. 
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2 Laplace-Beltrami operator and 2-dimensional hydrodynamics 
2.1 Complex and symplectic geometry 

Hereafter S will always denote a two dimensional compact boundaryless orientable manifold (hereafter called closed), 
endowed with a riemannian metric g = {,), and D C S a Jordan domain (a region with compact closure on a 
whose boundary consists of a finite number of connected components each diffeomorphic to a circle). The underlying 
Riemann surface structure for S is given by the atlas A of Laplace-Beltrami isothermic coordinates. Complex geometry 
keeps the notion of angle between tangent vectors u, v (thus the rotation operator J by 90 degrees of tangent vectors 
is well defined), but neglect lenghts |u| and areas Q(u,v). Symplectic geometry, on the other hand, keeps the area 
form but neglects the angles. Merging different symplectic forms with the same complex geometry produces conformal 
riemannian metrics. In 2-dimensions all metrics are automatically Kahler. 

Functions (usually velocity potentials or stream functions) produce vectorfields under the gradient and symplectic 
gradient operators grad (usually denoted V), sgrad := J o grad (also denoted V x ), 

#(•) = (grad^, •} = Q(sgrad(0), •) . (2.1) 

There is a "baby Hodge theory" in two dimensions: one identifies functions with two-forms via / <-> fQ and vectorfields 
with 1-forms via v «-> w = («,•)■ Let C be a closed curve enclosing a domain homeomorphic to a disk D C S. Orient 
C such that the frame n, t is positive, where t is the unit tangent vector and n the exterior normal. Let v be a 
vectorfield. The familiar formula 



(v,n)ds= / / div(u)fi (2.2) 

I C J J D 

makes perfect sense. The coordinate free definition of the divergence operator, sending vectorfields to functions, can 
be given as follows (applying Stokes' theorem and looking at the integrands): 

d{J(v), •) = div(u)0 . (2.3) 

Incompressible hydrodynamics means: div(u) = except for sources or sinks. On a compact surface S their divergences 
must add to zero. In potential flow (and always incompressible), v = — grad</!> for a velocity potential <f> which in general 
is not globally defined (it is a multivalued function in the Rieman surface sense). 

Define the Laplace operator by A = div o grad. It is a negative definite self-adjoint operator with respect to 
(/,<?) = J fgQ,. Incompressibility implies A(f> = (up to singularities). Thus the velocity potential <p(s),s £ S is 
harmonic up to singularities, and complex function theory of S (viewed as a Riemann surface) comes to fore. The 
conjugate harmonic function ip is called the stream function. Combined they form the velocity potential F — (f> + iip 
(in planar hydrodynamics, u + iv = — F'(z)). Clearly v is tangent to the level lines of ip, and can be described 
symplectically from ip via 

dtp = Q(«, •) <4> v = sgrad(?/>) = -grad(0) . (2.4) 
The circulation § c (v,t)ds of v around a closed curve C = dD is an object of fundamental importance. Using Stokes' 
theorem it can be transformed into a double integral, 

(v,t)ds= <j> (J(v),J(t))ds= (L (J(v),-n)ds = - / / div(Jv)n 
c Jc Jc J Jd 

But J{v) = J(Jgradi/>) = — grad^- Hence 

j> (v,t)ds =11 A^n (2.5) 

J C J J D 

A coordinate free definition of the surface rotational or vorticity ui = rotv is then 

d(v, •) = rot(v)fi . (2.6) 

In passing, we mention that discrete versions of all these objects have been proved useful in computer graphics 
and flow visualization, see e.g. [23], [53] . 
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2.2 Poisson's equation and Green functions for the Laplace operator 

The vorticity distribution ui forms "the sinews and muscles of fluid motion" (Kuchemann, [33], 1965). No wonder 

2-dimensional hydrodynamics is governed by Poisson's equation (inverting the Laplace operator). Perhaps the noblest 
equation of them all, it is the holy grail of every mathematician, pure or applied: 

Atjj = uj , ui G some class of functions . (2-7) 

On a closed surface the source term must average to zero: 

un = . (2.8) 



s 

[ Proof. J J s uiQ = J J s AipdS = J J s V • VipdS = <f gs (Vi/j, n)ds — 0, since there is no boundary. Alternatively, take 
a small curve C enclosing a disk D. The circulation of the velocity field can be computed as the double integral of 
the vorticity on D and minus that integral on S — D.] 

In two dimensions vorticity is a "material property", L v uo = 0, and this translates into 

^ + (v, gradw) = , v = sgrad?/) . (2.9) 

Choosing an appropriate class of vortex functions is an important mathematical problem. In many concrete 
examples, the "vorticity ountour map" has sharply defined islands of high vorticity (positive or negative) in a sea of 
roughly constant (often zero) vorticity. In time, such regions can blend or develop filaments. The evolution of the 
system (|2.7I I2.9p may reach a stand-off by the appearance of singularities. 

At any rate, given Kj , 1 < j • < N, take geodesies disks Dj of equal but very small area Area(Dj) = e 2 centered 
at Sj and define 

u(s) = % , s £ Dj (2.10) 

and the suitable constant elsewhere. One can prove that an initial state consisting of small nearly circular vortex 
cores remains stable. 

A vortex pair is the case N = 2 with opposite vorticities. Their ubiquous presence in fluid flows is fascinating. In 
section [4731 we give a simple proof of a conjecture by Kimura [35]: a dipole describes a geodesic on S (zooming across 
the fluid very quicky). 

A monovortex consists on a single vorticity island on a sea of distributed counter vorticity (uniformly, say) so that 
the average vorticity is zero. By linearity, any system of vortex patches can be decomposed on a weighted (positive 
and negative) sum of monovortices. The following gedanken experiment may help. Consider a fluid at rest occupying 
the closed surface S. All of a sudden, stir around a given point s , producing a vortex singularity there. What 
would be the resulting flow? Since there is no distinguished point s, to be taken as a counter-vortex, the opposite 
vorticity should be distributed as homogeneously as possible on the surface. [ Vorticity is an area related concept (via 
Hodge theory, vorticity is thought as a two-form), and hence the area element of S should be used to perform this 
homogenization.] Note that a distinguished s„ may exist if S has symmetries, for instance if s is the south pole of a 
surface of revolution; then s* is the north-pole; those situations would be quite unstable, thouglQ. 

Point vortex systems are meant as finite dimensional ODE aproximations for Euler's equations (geodesic flow on 
the group of area preserving transformations of S) or the equivalent PDEs (|2.7I2.9[) . assuming concentrated vorticities 
at points Sj(t), 1 < j < N. As e — * 0, how to obtain a suitable ODE approximation for the PDE system (|2.71 12.9)) ? 
This is done in two steps. Firstly, one needs the stream function for a marker particle on a flow generated by bound 
point vortices. The second step is a desingularization procedure in order to remove the influence of each vortex on 
itself. Then we let them "dance" under each other's influence and its own drift. 

The first step is achieved with Green functions od the Laplace-Beltrami operator. 



2 The lack of unicity in some Euler flows is a known possibility; Joe Keller's "tea pot effect" is a sometimes dramatic 
example [32] . 
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3 Point vortex equations 

Definition 1 G — G(s,g) * s the Green function of A g , 

A~ 1 u>(s) = / G(s,r)u{r)dA(r) , A -1 Aw = lu - \ / ujdA . (3.1) 

Jm a Js 

G(s,s ) is a smooth function outside the diagonal, characterized by the following properties ]50f : 

AG(s, s ) = - Are ^g^ ( s / s °) > G(s,s D ) - log d(s,s )/27r bounded, J G(p, q)Q{q) — , G(s,s a ) — G(s a ,s). 

(3.2) 

Symmetry is to be expected from this formal calculation: 

G ao (s)= / / G So ( Sl )AG s ( Sl )dsi = I I VG ao -VG s - / G ao d n G s = [ [ VG So ■ VG a . 

J JS J J J dS J J 

Hydrodynamically, ip — G(s,s ) fulfills all the requisites to be the stream function of a unit monovortex at s . 
In the electrostatic interpretation, G(p; q) is the potential at p of a positive unit point charge at q, on a conductor 
with total charge zero (the background negative charge is uniformly distributed). A test particle s with positive but 
negligible charge on the field generated by a unit charge will move according to 

s = gradG(s; s D ) , s / s D . (3.3) 

while a fluid particle (marker) s on the flow generated by an unit strength bound vortex s will move according to 

s = sgradG(s;s ) , s ^ s . (3.4) 

3.1 Eppur si muove: vortex drift under its own stream 

Newton's third law precludes the action of the material particle at s on itself, but the rest of the universe S — s con- 
spires to impinge a "collective reaction" on s D and makes it drift. In order to capture this motion, a desingularization 
procedure is needed. To make that long story short, a single vortex will move according to 

Prescription 1 

s = sgradi?(s D ) , R(s„) = lim G(s,s g ) - logrf ( s ' s °) (3 5) 

where d is the geodesic distance between points as measured with g. 

The explanation for this rule is given in Proposition [1] below. We were filled with a religious feeling by the fact that 
the regular part 7? is a notable object from the geometric function theory of the Laplace-Beltrami operator, known 
as Robin's function of g ([51], [H]). Vortex drift belongs to mainstream mathematics! 

Remark 1 Jean Steiner 159/ has shown that for any metric g on S 2 

R g (s) - i- (A 9 - 1 A- 9 )(s) = -^traceA; 1 (3.6) 

where K is the Gaussian curvature. For higher genus the RHS is no longer a constant and its fiutuation does not 
have a simple geometric interpretation. 



Remark 2 Marker's equation [3.$ becomes a 1 1/2 degrees of freedom system, by inserting s (t), a solution of f3. 51) 
Chaotic marker equations are possible even for one vortex systems. 
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A basic observation in this paper is that the core energy argument in [27] geometrizes nicely: 

Proposition 1 The desingularization rule 13. 5j) encodes the core energy renormalization argument. 

Proof. Consider a small Gauss system of normal coordinates around s . The gradient gradd(s, s ) is a unit vector 
along the geodesic rays. The symplectic gradient sgrad d(s, s ) is obtained by composition with J and is therefore 
tangent to the geodesic circles. Its flow leaves the geodesic disks invariant, and the same is true for any function of 
the distance, in particular losing = logd(s,s ). For log d(s,s ) the vectors rotate with speed inversely proportional to 
d(s,s ) on a geodesic circle. A simple back of the envelope reasoning shows that the kinetic energy confined in this 
blob is logarithmically infinite. It can be removed (in physics jargon, renormalized) provided only a small quantity 
of kinetic energy crosses the boundary of a geodesic disk in finite time. We invoke |27| : energy diffusion "can be 
neglected in the limit as the radius of the ball tends to zero" , section 5). This is true precisely due to the the fact 
that behavior of Greens's function G(s,s ) — d(s, s )/2n is bounded. 



3.2 Collective vortex motions: Hamiltonian description 

Collective motions are obtained under 

Prescription 2 Each vortex moves on the stream generated by all the others together with its own drift term. 

Rewriting this prescription in Hamiltonian form is straightforward. In view of the symmetry properties of Green's 
functions, C.C.Lin's [45] arguments geometrize directly. 

Theorem 1 The dynamics of N point vortices Sj of strenghts Kj moving on S is a Hamiltonian system on S x ... x S 
with symplectic form 

^collective ~ ^ ^ ^J^(^j) (^ - 7) 
l<j<JV 



and Hamiltonian 

N 



H = }_^ KiKjGisi, Sj ) + -k* R{ Si ) . (3.8) 

i<Cj i— 1 

where the self interaction terms are given by Robin's function 

fl(,) = lim G(r,,)- lQg(d o (r ' s)) . (3.9) 

Remark 3 

i) The key aspect for deriving vortex equations is the geometric desingularization prescription f3.9V . As we just 
observed, it can be solidly justified via the core energy method, see section 5 of Flucher and Gustafsson ( j27j , 1997). 

ii) We omitted the possibility of collisions. See 129$ . \30$ for collisions and regularizations for 3 and 4 vortex motions 
on the plane. 

There is an immediate extension to vortices with mass that may be relevant due to current interest in Bose- Einstein 
condensates, see [40], [13] . [54]). These systems will exhibit slow-fast Hamiltonian phenomena [48] . 

Corollary 1 Let each Sj have a mass irij besides its vorticity Kj, contributing with a kinetic energy (p 3 j,Psj)> 
where the bracket denotes the the induced inner product in T* S via the Legendre transform of g. The dynamics in 
T* (S x .... x S) is described by the Hamiltonian system 

H = ^ (p Sj ,p Sj ) + ^ KiKj G(Si,Sj) , ^collective = ^can + ^ Kjtt(Sj) (3.10) 

where Q. ca n = "2 ^Pi ^ " 15 ^ e canonical 2-form ofT*(S x ... x S). 
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We now address the following question. Let g — h 2 g a conformal metric. How do the vortex problems on (S, g) 
and (S, g) compare? Of course, the symplectic form deforms as Q, — » h 2 Q in each coordinate. Can we find the new 
Hamiltonian in terms of the old one? We need transformation formulas both the Green and Robin functions. 

3.3 How Green and Robin functions change under conformal transformations 

Let g — h 2 g a conformal change of metric. Green's function Gg(s,s ) for the Laplace- Beltrami operator on a closed 
surface is not conformally invariant. It must change, because a background uniform vorticity is an area dependent 
notion, except when it is zero. In fact, if the background vorticity is zero, the following "naive conformal rule" holds: 



Proposition 2 In a zero vorticity background, in order to find the g-regularized stream function for a vortex, subtract 
log(fa)/27r from its g-regularized stream function, g = h 2 g. 

Proof. 

V'reg (So) = lim ^conformal (s) — log rf(s, S )/27T . 

a — >-s 

/ d(s,So)\ 

log(d(s, s )) = log d(s, s ) ' ° ~ log(d(s, s )) +logft(s ) for s near s . 



d(s,s ) 

However, if, and this is a big if the sum of vorticities does not vanish, there is a constant nonzero background 
counteracting vorticity in S — {si, ...,s n }. Since what is constant is a metric dependent notion, this complicates the 
conformal transformation rules. Using Green's function for the Laplace-Beltrami operator on a closed surface becames 
unavoidable. The correct transformation formulas are as follows 



Lemma 1 (|50j) 

G(s,s )-G(s,s ) = -\(A- 1 h 2 (s) + A- 1 h 2 (s )) + -^ [ h 2 A^h 2 n (3.11) 
A A J s 

R(s) = R(s)- i/i- jAg 1 h 2 (s) + i J h 2 A^h 2 Q, (3.12) 

A beautiful theory has been recently developed around Robins's function and spectral invariants of conformal classes 
of metrics, see |59| . For our purposes the last term in both formulas can be dropped, as they just add constants. After 
some algebra juggling, one gets 

Theorem 2 Under a conformal change of metric g(s) = h 2 g(s), the symplectic form deforms as f2 — > h 2 Sl in each 
coordinate, and the Hamiltonian becomes 

^ JV JV JV 

H = H — —} K 2 ~iog(h(sA) > KtA~ 1 h 2 (sA , «=> m . (3.13) 

4-7T * — ' A £ — ' £ — ' 

e=i i=i i=i 

Remark 4 A map F : (S,g) — > (S,g) is conformal when 

g{dF s -v s , dF s ■ w s ) = h 2 (s)g(v a ,w s ) . (3-14) 

The left hand side induces a metric gF(v s ,w s ) in S, hence g~F = h 2 g and Theorem^ applies. The vortex dynamics 
(Q,,H) on S x ... x S corresponds via F to the dynamics of vortices in S x .... x S. 
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For completeness, we present the Proof of lemma[T] following [51] . The trick is average over Cl twice using the "axioms" 
in Definition [T] Consider the functions 

G(s, s ) - G(s, si) = E(s; s , si) , G(s, s a ) - G(s, si) = E(s; s a , si) . 

Both E and E are harmonic up to (+) and (-) log singularities at s and si so they differ by a constant c. To find 
this constant we do the fl(s) = h 2 (s)Q(s) average of 

E — E = {G{s, So) - G(s, si)) - (G{s, s a ) - G(s, si)) = c . 

The first two terms drop out while the last two give 

A _ V(*i) - A~V( So ) = cA . 

Hence 

(G{s, so) - G(s, si)) - (G{s, so) - G{s, s x )) = A" 1 ^^)^ - A~ 1 h 2 (s )/A . 
Again, average this expression, but now over f2(si). We get 

G{s,So)A-0-G{s,s o )A + A- 1 h 2 (s) = / A" 1 /i 2 (s 1 )/i 2 (si)n(s 1 )/A - A _1 /i 2 (s ) . 

Js 

This gives the transformation formulas G — > G. The transformation formula R — > _R follows by the limit s — > s in 

g (s , So) - Mfeil) = ( G(S] So) _ G(S; So)) + f G(S] So) _ iQg^.^) ^ _ iog(d( a , ao )M 8 , ao )) 
3.4 Jordan domains, Schottky doubles 

The "naive rule" of Proposition [2] works in the case of a Jordan domains D C S, regardless of the vorticity sum. The 
basic ingredients for vortex motion are the hydrodynamical Green functions: 

Definition 2 A hydrodynamical Green function Gd(s;s ) is a real harmonic function of p G D — s , with (-) loga- 
rithmical singularity at s — s , symmetric in s,s , continuous in D — s . G is constant on the boundaries and with 
desired circulations around them (topology imposes relations). 

C.C. Lin's result on the added term in the Hamiltonian under a conformal transformation is easily interpreted in 
terms of Lemma [2] the "naive rule". Let / : (D, S,g) — > (D, S,g) a conformal map between two Jordan domains, 

g(df(s) ■ u,df(s) ■ v) = h 2 (s)g(u,v) . (3.15) 



Proposition 3 Let Gd{s; s ) be the hydrodynamical Green function for D . Then: 

Gp(s,s~o) = Gd{s;s ) (3.16) 
The regularized Green function receives the naive correction, 

9(so)=g(so)- 1 -^ . (3.17) 
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Note that the sign in equation (4.6) of [28] must be corrected from + to -. 

For domains conformal to a planar region Green functions are classified into two types. First type, more employed 
in electrostatics, means that G — on all boundaries, possibly with different circulations around each of them. Second 
kind, or modified Green functions are more commonly used in hydrodynamics. One prescribes zero circulation around 
the boundaries, except at one of them; the values of Gd = const, on the boundaries may differ, one of them being 
normalized to 0. The existence of such Green functions For multiply connected planar domains is due to Koebe, and 
was used by C.C.Lin in his seminal vortex paper [45]. 

Specially exciting is the fact that recently Crowdy and Marshall ([18].|19j.|21], [22] ) brought to "implementational 
fruition" the task if obtaining Green functions in canonical domains. We believe hydrodynamical Green functions 
should exist on arbitrary Jordan domains on Riemann surfaces - the difference from planar regions is that in general 
the domain will have "handles" (like a juice jar). 

Remark 5 An interesting viewpoint is to consider the Schottky double Sd — D + D of a Jordan domain D, S is 
called a symmetric Riemann surface. Consider a metric on S such that the reflection is an isometry. If at initial 
time each vortex on D corresponds to a symmetric one in D, with opposite vorticity, then the corresponding vortices 
remain symmetric for all time. N-vortex motion on D is the restriction of the corresponding 2N vortex motion on S 
(using the corresponding Hamiltonian from I). Thus Schottky double construction is reminiscent of Thomson's image 
vortex method. Inspite of the fact that a given metric in D in general will not extend to an reflection-isometric one, 
since the total vorticity is zero, only the Riemann surface structure is needed. In fact, Schottky-Klem prime functions 
have been proved useful to finding hydrodynamical Green functions on multiply connected planar domains ('201). 



Recall that when the sum of vorticities does not vanish, using Green's function for the Laplace-Beltrami operator on a 
closed surface is unavoidable, because there is a constant nonzero background counteracting vorticity in S— {si, s n }, 
and this complicates the conformal transformation rules. 

4.1 Riemann surface Green functions 

The case where ^2 K i — h as a special feature: the collective vortex stream function is the same for all conformal 
metrics. Of course, this is the basis of conformal mappping methods in applied hydrodynamics. Mathematically, this 
follows from the special two-dimensional property 



Hence, irrespective of the chosen metric in its conformal class, all Laplace-Beltrami operators annihilate the same 
functions, logarithmic singularities allowed. We need only the complex structure A: residues at poles are conformal 
invariants, etc. 

Lset Gs = G(s; s , si) be the two point Green function for (S,A), whose existence for closed Riemann surfaces 
S can be used as the starting poin10 for the theory of meromorphic functions and differentials and uniformization 
theorems, see e.g. Weyl ([62] ,11. 13, potential arising from a doublet source) or [58]. G is the unique (up to a constant) 
real harmonic function for s £ S, except for a (+) logarithmical singularity at s and a (-) logarithmical singularity 
at si. 

53 =i K i = implies that the collective vortex stream function can be given in a conformally invariant way: 



4 When the total vorticity vanishes 



g — h 2 g =>• Ag = h 2 A B . 



N 




(4.1) 



3 



It is possible that Riemann [55] and Klein [3S] insights came from attended Helmoltz lectures. 
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where s* is any chosen point in S. Irrespective of the metric, the background vorticity of every term in zero except at 
the poles Sj, s* . Moreover, the choice of s* is irrelevant as the infinite vorticities there cancel out because ft, = 0. 

[Unfortunatelly this is not always a free lunch, since there are some drawbacks. One is that Riemann surface 
two point Green functions have no symmetries in its arguments, which makes it impossible to use G(s; Sj , s*) for a 
collective Hamiltonian formulation. Secondly, there is no direct interpretation for the desingularization procedure in 
terms of invariant objects such as Robin's function. One can also add a practical objection, that it seems to be equally 
hard finding (in practice) Gs, g {s,s ) or Gs(s; s , si). In fact, Gs(s; s D , Si) = Gs, g {s,s ) — Gs, g (s,si) + const, where 
in the right hand side any metric can be used.] 

At any rate, for any vortex problem in which the background vorticity is zero, the "naive conformal rule" (Propo- 
sition [2)| holds, and this fact may simplify matters considerably. 



4.2 Surfaces conformal to S 2 

Take the unit sphere S 2 C 5ft 3 as the "concrete model" for S (thought of as an abstract Riemann surface). Any metric 
can be written as g — h 2 g a , where g is the canonical metric on the sphere. Stereographic projection F : S 2 — > C 
from the north pole s* = (0, 0, 1) to the equatorial z-plane yields [35] 

g = h 2 \dz\ 2 , ho = 2 . (4.2) 
1 + \z\* 

We can use the z-plane to study the dynamics. If a particle passes through the north pole, in C it will disappear 
in the infinite but reappear instantaneously from another direction, but ee may assume that this is not happening 
during the time frame we are working witlQ. In z-plane, the Green function is log(«)/27r. Let h ■ h the combined 
conformal factor from (S 2 ,g) to C. Because of zero total vorticity, vortex dynamics on (S 2 ,g) can be studied in C 
with the symplectic form 

ft = ^2 Kjh(xj,yj) 2 h (xj,yj) 2 dxj Ady 3 , (4.3) 

3 

and Hamiltonian given by Theorem 13.131 (we will omit the factor 2n): 

1 ' ^1 
H= -^Kj/t^logd^ - zi\ 2 ) ~2J -n 2 m \og{h {x m ,y m )h{x m ,y m )) . (4.4) 

3,1 P=l 

Here ' means that the j = £ terms are omitted in the summation. Now we rewrite each re„ as — K m (^2 p v ^ m k v 
expand and regroup terms. We also use the following interesting formula: 



Lemma 2 (|16j) 

si - s 2 | 2 = ho(si)h (s2)\zi - z 2 \ 2 . (4.5) 
Here the distances are the euclidian distances in 5ft 2 and 5ft 3 . 

Pulling back to the sphere, we get: 

Proposition 4 For a system of N point vortices in the sphere S 2 with a metric g = h 2 g (g a is the standard metric), 
and such that the total vorticity vanishes, the dynamics is governed by 



Kjh 2 (si)iu(si) , H = i ^ KjKe log (h(sj)h(st)\s 3 - - s e | 2 ) (4.6) 

3 3,1 

where \sj — se\ is the euclidian distance and ui is the area form of the sphere. 



4 Esse est percipi, stealing the motto of Bishop Berkeley, ironically the greatest foe of Calculus lovers. 
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4.3 Proof of Kimura's conjecture on dipole motion 

The Hamiltonian for a vortex pair with opposite vorticities can be written as 



H = -k i + k B(si,S2), B(si,s 2 ) = 

Z7T 



R(si)+R(s 2 ) 



- G(s u s 2 )- 



logd(si, S2) 
2^ 



(4.7) 



Let k = O(e) and initial conditions d(si(0), $2(0)) = O(e). Kimura's conjectured that as e — > the vorticity center 
moves on a geodesic. The truncated system taking only the first term in (14. 7[) yields a system of ODEs on D(s ) xD(s ), 
the product of two copies of a geodesic ball centered on s 6 S, given by 



li = -/tsgrad ai logd(si,s 2 ) , s 2 = KSgrad a2 logd(si,s 2 ) 



(4.8) 



We claim that this equation represents the dominant O(l) geodesic dynamics plus an 0(e 2 ) perturbation. Due to its 
symmetric form, B(si,s 2 ) also produces an 0(e 2 ) perturbation. We take Gauss coordinates ([60], [23] ) 



ds 2 = du + G(u, v)dv 2 , G(0, v) = 1 , ^- G(u, v) = 

k»=o 



(4.9) 



The u-curves form a field of geodesies, and the central i>-curve (for u — 0) is also a geodesic. Suppose s corresponds 
to some v and u = 0. In these coordinates, si(0) = (— e, u) and S2(0) = (e, v). Clearly, at t — 0, si,2 will be tangent 
to ti-curves, and in the Gauss coordinates, 



t)i(0) = K/2e- 



•u 2 (0) = K/2e 



1 



mi(0) = u 2 (0) = 



In the limit as e — * (|4.9p . with k = 2ft e, k = 1, we get 

wi(0)=wa(0) =l + 0(c) , 6i(0),ua(0)=O(e) . 
Note that G(w, u) = implies that (|4.8[) does not have an O(e) term. This concludes the proof. 



(4.10) 



(4.11) 



5 Final Remarks 

5.1 Robin's function and vortex drift 

A single vortex s = sgrad_R(s D ) describes contour lines of Robin's function, so it would be interesting to understand 
its generic properties. For metrics on the sphere R satisfies (see also (|3.6I 59 ) 

AR = h 2 (K E -4i:/A(E)) (5.1) 

where A is the standard Laplacian on the sphere. 



5.2 Batman's function and integrable vortex pair problems on Liouville surfaces 

If S embedds in R 3 having an axis of symmetry n, the momentum map of the S 1 action on S x ... x S with 
SI = ^Kjf2 s . is J — Q^KjSj) ■ n. In particular, a vortex pair on a surface of revolution is integrable, with second 
integral / = (si — s 2 ) ■ n. Can we find nontrivial examples of integrable vortex pairs? As it is well known, Jacobi 
has shown that the geodesic flow on the ellipsoid is integrable. In view of Kimura's theorem, vortex pairs extends 
geodesic problems, so the following question is natural: 

Question 1 Is the vortex pair problem on the triaxial ellipsoid (with opposite vorticities) integrable? If it is not, can 
one prove it analytically by Melnikov's method applied to the geodesies passing by the umbilics? 
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Conformal maps from the triaxial ellipsoid to the standard sphere are known ([17]! HZL [SS])- The conformal factor 
h 2 can be retrieved from these papers, so Proposition U can be implemented. Alternatively, using sphero-conical 
coordinates (see, eg. [39]) allows obtaining h 2 from scratch with two elliptic integrals. More generally, Liouville 
surfaces (see e.g., [9]) are those for which the geodesic flow is integrable. For which Liouville surfaces is the vortex 
pair system integrable? In this programme we have found some preliminary results [?T] for (|4.7|l showing how the 
geodesic system is perturbed for d(si,S2) = 0{e). The 0(e 2 ) perturbation term is 



which could be called, for the sake of justice, Batman's function. 

5.3 Time dependent problems 

Either on closed surfaces or Jordan domains it is natural to consider time dependent problems, subject to the natural 
constraint (in view of imcompressibility) that the total area with respect to the metrics g(t) does not change. For the 
sphere and disks Riemann mapping theorem guarantees that the change is always conformal. In the general case the 
complex structure may also change, an interesting complication that should be addressed. 

5.4 Higher dimensional extensions 

Alan Weinstein called our attention that the statement of Theorem [1] makes sense on a compact Kahler manifold, 
replacing the two dimensional objects by their higher dimensional analogues. In particular, Calabi-Yau manifolds 
are very fashionable objects nowadays. Would the generalization of (|3.71 13 . 8f) represent singular solutions to Euler's 
equation for the hydrodynamics on (S,g), i.e. the geodesic flow on SDiff(5'), in the spirit of [3] and |46|? Or will that 
construction be more in line with field theorists approaches, see e.g., Baez [3]? 

5.5 Vortices for a proof of Riemann hypothesis? 

One of the approaches for a proof of Riemann hypothesis is searching a connection between the notrivial zeros of 
Riemann's zeta function with some quantum mechanical problem (Berry and Keating, [5])- The quantum version of 
geodesic motion on (S,g), i.e., the free particle, leads to the spectrum of minus the Laplace-Beltrami operator. As a 
vortex dipole moves on geodesies, why not look for a proof quantizing a vortex problem, which is asymmetric with 
respect to time reversal? 

"Sir Michael Berry proposes that there exists a classical dynamical system, asymmetric with respect to time 
reversal, the lengths of whose periodic orbits correspond to the rational primes, and whose quantum-mechanical 
analog has a Hamiltonian with zeros equal to the imaginary parts of the nontrivial zeros of the zeta function. The 
search for such a dynamical system is one approach to proving the Riemann hypothesis (Daniel Bump, [14])." 

We conclude this note on a highly speculative tone: 

Conjecture 1 A connection may exist between the zeros of Riemann's zeta function and the quantization of a 3/2 
degrees of freedom vortex monopole problem on some compact Riemann surface with metric g, 




s — sgrad G(s,s ) , s a — sgrad_R(s Q ) , 



or with a vortex pair 



si = - 



Ksgva.d Bl H(s 1 ,s 2 ) s'i = /tsgrad S2 .ff(si, s 2 ) , H 



= — K 



2 logd(si,s 2 ) 
2tt 



+ k 2 B(si, s 2 ) . 
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